Some Global Characteristics of the Galactic Globular Cluster System 



Sidney van den Bergh 

Dominion Astrophysical Observatory, Herzberg Institute of Astrophysics, 
National Research Council of Canada, 5071 W. Saanich Rd., Victoria, British Columbia, 

V9E 2E7, Canada 

Sidney . vandenbergh@nrc-cnrc . gc . ca 
ABSTRACT 

The relations between the luminosities My, the metallicities [Fe/H], the Galacto- 
centric radii R, and the central concentration indices c of Galactic globular clusters are 
discussed. It is found that the most luminous clusters rarely have collapsed cores. The 
reason for this might be that the core collapse time scales for such populous clusters 
are greater than the age of the Galaxy. Among those clusters, for which the structure 
has not been modified by core collapse, there is a correlation between central concen- 
tration and integrated luminosity, in the sense that the most luminous clusters have the 
strongest central concentration. The outermost region of the Galaxy with R > 10 kpc 
was apparently not able to form metal-rich ([Fe/H] > —1.0) globular clusters, whereas 
such clusters (of which Ter 7 is the prototype) were able to form in some nearby dwarf 
spheroidal galaxies. It is not yet clear how the popular hypothesis that globular clus- 
ters were initially formed with a single power law mass spectrum can be reconciled with 
the observation that both (1) Galactic globular clusters with R > 80 kpc, and (2) the 
globulars associated with the Sagittarius dwarf, appear to have bi-modal luminosity 
functions. 

Subject headings: globular clusters: general 



1. INTRODUCTION 

Globular clusters are among the oldest objects to have formed in the Galaxy. They therefore 
provide valuable information on the early evolutionary history of the Milky Way System. The 
present paper discuses the interrelationships between integrated cluster magnitude My, cluster 
metallicity [Fe/H], Galactocentric distance R, and the central concentration index c = log(rj/r c )(King 
1962). The basic observational data, which are listed in Table 1, were drawn from Harris (1996), but 
updated from new information that is compiled at http://physun.physics.mcmaster.ca/Globular.html. 
A discussion, in some ways similar to the present one (but based on older and less complete data), 
was previously published by Djorgovski & Meylan (1994). That paper also gives extensive references 
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to earlier work on this subject. Traditionally (Kinman 1959; Zinn 1985) the Galactic globular 
cluster system is regarded as consisting of a metal-rich disk and a metal-poor halo. However, more 
recently (Minniti 1995; Barbuy et al. 1999; Cote 1999; van den Bergh 2000a) it has become clear 
that the metal-rich clusters actually constitute a bulge, rather than a disk, population. Among the 
questions to which answers will be sought in the present investigation are the following: (1) Is the 
Galactic globular cluster system old enough for even the most massive globular clusters to have 
collapsed? (2) Does the structure of pre-collapse clusters depend on their mass or environment? (3) 
Do metal-poor globular clusters with [Fe/H] < — 1.0 exhibit a radial metallicity gradient over the 
range —0.5 < log-R(kpc) < 2.5? (4) Do metal-rich clusters with [Fe/H] > —1.0 show a radial abun- 
dance gradient over the range —0.5 < logi? < 1.0? (5) Did the metal-rich globular clusters in the 
outer reaches of the Galaxy have a different origin from the metal-rich globulars with R < 10 kpc 
that are embedded in the main body of the Galaxy and its inner halo? Questions (4) and (5) above 
may be rephrased by asking if the Galaxy contains physically distinct metal-rich and metal-poor 
cluster populations. The presence, or absence, of metallicity gradients within these two populations 
might throw some light on the evolutionary history of each of these two populations. Finally it is 
noted that attempts to answer some of these questions are intrinsically uncertain because of the 
relatively small number of globular clusters associated with the Galaxy. 



2. CORRELATIONS 

All six possible permutations of the relations between My, [Fe/H], R and c are plotted in 
Figure 1. A brief discussion of observational data on these relations is given below. 



2.1. [Fe/H] versuslogR 

This plot shows the well know dichotomy between the metal-poor clusters with < [Fe/H] 
—1.5 and the metal-rich clusters with < [Fe/H] > ~ —0.5. All but three of the metal-rich 
{[Fe/H] > —1.0) globulars are located at R < 10 kpc, i.e., almost all metal-rich clusters are located 
in the inner region of the Galaxy. Among these exceptions two (Terzan 7 and Palomar 12) seem to 
be associated with the tidally disintegrating Sagittarius dwarf galaxy (Ibata et al. 1994; Irwin 1994). 
The third exception is the cluster Palomar 1, which appears to be a "young" globular (Richer et al. 
1996) that might have formed in a no longer extant dwarf spheroidal (van den Bergh 2000b). In 
summary, 43 out of 46 (93%) of metal-rich Galactic globular clusters with [Fe/H] > —1.0 are old 
objects located at R < 10 kpc. The three remaining clusters are thought to be younger objects that 
formed (or may have formed) in dwarf spheroidal galaxies. Is there a radial metallicity gradient 
among metal-poor globular clusters? Such a gradient would be expected (Eggen et al. 1962) if the 
young Galactic halo was able to enrich itself in heavy elements while it was contracting. On the 
other hand, such a gradient would not be expected for a more chaotic assembly of the halo, such 
as has been proposed by Searle & Zinn (1978). Data on the radial distributions of moderately 



-3- 



metal-poor (-1.60 < [Fe/H] < -1.00) and very metal-poor (-2.40 < [Fe/H] < -1.60) clusters 
are collected in Table 2. Separate information is provided for (1) all metal-poor clusters, and (2) 
with probable companions of the Sagittarius dwarf excluded. Inspection of these data hints at a 
possible metallicity gradient in the expected sense, i.e., with a small excess of the highest metal- 
licity subgroup at small Galactocentric radii. However, a Kolmogorov-Smirnov test shows only an 
88% probability that this effect is real for the entire sample, and an 82% probability that it is real 
for the sample from which probable companions to the Sagittarius galaxy have been excluded. In 
summary, it is concluded that the population of metal-poor Galactic halo clusters is too small to 
establish with certainty if they exhibit a metallicity gradient. For metal-rich globular clusters with 
[Fe/H] > — 1.0 there is no statistically significant evidence for a correlation between [Fe/H] and 
R. Finally, it is noted that we presently do not understand how some nearby dwarf spheroidal 
galaxies were able to form quite metal-rich globular ([Fe/H] < —1.0) clusters. This contrasts with 
the situation in the outer (R > 10 kpc) regions of the Galaxy that were apparently never able to 
form similarly metal-rich globular clusters. This result suggests that the specific globular cluster 
frequency in the Galaxy may have exhibited a steep decline with increasing metallicity, whereas 
the specific cluster frequency may have been less sensitive to metallicity in (some) dwarf spheroidal 
galaxies. 

2.2. [Fe/H] versus My 

The metallicity distribution of Galactic globular clusters does not appear to depend on luminos- 
ity. Kolmogorov-Smirnov tests show no significant differences between the metallicity distributions 
of objects brighter and fainter than My = —5.0. The same conclusion holds if the sample is divided 
at My = —6.0. This result is somewhat surprising because one might have expected faint metal- 
rich clusters (which mainly occur in high density regions) to have been preferentially destroyed by 
bulge and disk shocks, or eroded by tidal stripping. On the other hand the most massive (luminous) 
clusters are expected to be survivors e.g., Gnedin (2002). 

2.3. [Fe/H] versus c 

The Figure shows little (or no) evidence for a dependence of metallicity on the central concen- 
tration of index c (King 1962) of Galactic globular clusters. In particular the fraction of globulars 
with collapsed cores (c « 2.5) does not appear to depend strongly on metallicity. 

2.4. c versus log R 

The overwhelming majority of clusters with collapsed cores, i.e., those with c ~ 2.5, are found to 
be located at R < 10 kpc. The reason for this is, no doubt, that globular cluster half-light radii tend 
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to decrease with decreasing Galactocentric distance (van den Bergh & Morbey 1984). As a result 
the typical relaxation times (and core collapse time scales) of clusters near the Galactic center are 
significantly shorter than they are for globulars at larger Galactocentric distances. Among metal- 
rich {[Fe/H] > —1.0) clusters the central concentration index c correlates with Galactocentric 
distance in the sense that the most metal-rich clusters have the strongest central concentration of 
light. For metal-rich clusters with c < 1.00, < logi? > = 0.69 ± 0.07, which is significantly larger 
than < logi? > = 0.34 ± 0.10 for the metal-rich clusters with c > 2.00. A similar correlation for 
metal-poor clusters appears to exhibit more scatter than does that among the metal-rich Galactic 
globulars. The observation that clusters with low central concentration, on average, have larger 
Galactocentric distances than do those at larger distances is (see Section 2.5) due to the fact that 
most faint clusters with My > —6.0 have low central concentrations (c < 1.0), whereas the majority 
of clusters 1.5 < c < 2.4 are more luminous than My = —6.0. 



2.5. c versus My 



The structure of clusters may be characterized by their central concentration index c. A plot 
of c versus My exhibits a clear dichotomy between clusters that have collapsed cores with c ~ 2.5, 
and those with smaller c values that do not exhibit such collapsed cores. It is of interest to note 
that only one Galactic globular cluster with My < —8.0 is known to have a collapsed core. This is 
NGC 7078. The reason for the observed paucity of luminous clusters with collapsed cores (c = 2.5) 
might be that the core collapses time scale, which scales as N/lnN, is too long for the most massive 
clusters to experience core collapse within a Hubble time. However, a possible argument against 
this view is that few of the clusters for which Harris (1996) lists half-mass relaxtion times have 
relaxation times longer than 5 Gyr. A comparison between the frequency distributions of clusters 
with collapsed cores (c > 2.4), and clusters of intermediate central concentration (1.5 < c < 2.0), 
is shown in Table 3. A Kolmogorov-Smirnov test shows that there is a 99.1% probability that 
the intermediate concentration clusters, and the collapsed core clusters, were not drawn from the 
same parent luminosity distribution. A plot of c versus My, for those clusters that do not have 
collapsed cores, shows a close correlation between central concentration and luminosity. [The exact 
value of the correlation coefficient is sensitive to the choice of clusters that are omitted from the 
calculation because they are regarded as objects that are on their way to core collapse.] Among 
such clusters lacking collapsed cores the most concentrated objects have the highest luminosities. 
This correlation was previously also noted by Djorgovski & Meylan (1994). Clusters on the "main 
sequence" that lack collapsed cores have central concentrations that range from c « 0.6 at My s=s —5, 
to c « 1.8 at My ps —9. It seems likely that this observed correlation between luminosity and central 
concentration (of clusters that do not have collapsed cores) is due to initial conditions at the time of 
cluster formation. A few clusters, such as NGC6397 and NGC6717 lie between the "main sequence" 
and the sequence of objects with collapse cores at c « 2.5. Possibly such clusters are presently 
moving toward, but have not yet attained, core collapse. In view of the rich diversity of the relation 
between c and My, that is shown in Figure 1, it is remarkable (van den Bergh et al. 1991) that the 
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half-light radii of Galactic globular clusters turns out to be independent of cluster luminosity. It is 
noted in passing that most of the globulars associated with the Sagittarius system appear to lie on 
(or close to) the "main sequence" for Galactic globular clusters in the My versus c diagram. 

2.6. logi? versus My 

In the My versus logR plane centrally concentrated clusters with c > 2.0 are all found to lie 
below the line 

logi? = -0.43-0.21My (1) 

i.e., luminous compact clusters mainly occur at small values of logR, whereas faint extended 
clusters mostly lie at large Galactocentric radii. Inspection of the My versus log R diagram shows 
a clear distinction between the luminosity distributions of globular clusters in the main body of 
the Galaxy and those at large Galactocentric distances. Globular clusters with R < 80 kpc exhibit 
a near Gaussian 1 luminosity function that peaks at My s=s —7.5, while the outer halo clusters 
have either My s=s —5 or My s=s —10. Surprisingly (van den Bergh 2000a, p. 228) the luminosity 
function of the globular clusters associated with the Sagittarius dwarf system also has a bi-modal 
luminosity function that appears to resemble that of the clusters with R > 80 kpc in the outer halo 
of the Galaxy. [If one includes the clusters that may be associated with the Sagittarius, then the 
luminosity function might perhaps be better described as being weighted towards objects of very 
low luminosity (mass)]. Van den Bergh (2000a, p. 229) has used a Kolmogorov-Smirnov test to 
show that there is only a 4% probability that the luminosity distribution of Galactic globulars with 
R < 80 kpc were drawn from the same parent population as those as those with R > 80 kpc. It is 
widely believed (Fall & Rees 1977) that the presently observed globular clusters are the survivors 
of an initially much larger population. Within the main body of the Galaxy the present luminosity 
function of globular clusters has a near Gaussian shape (Abraham & van den Bergh 1995). It is 
thought that this is due to the dynamical "erosion" (Ostriker et al. 1972; Gnedin 2002) of what was 
initially a power law mass spectrum. However, it is difficult to see how such a process could could 
produce the apparently bi-modal luminosity distribution of (1) the Galactic globular clusters at 
R > 80 kpc, and (2) of the globular clusters that appear to be associated with the Sagittarius dwarf. 
Perhaps outer halo clusters such as N2419 (My = —9.6) and Pal.3 (My = —5.7) actually belong 
to separate populations with very different evolutionary histories. By the same token the cluster 
N6715 (=M 54) with My = —10.0, which is located at the center of the Sagittarius dwarf, may 
have had a different evolutionary history than that of the (and non-central) Sagittarius companions 
Ter 7, Ter 8, and Arp 2, all of which are faint and have My « —5. 



lr The "tail" of faint clusters with My > —4.0 may represent objects that were initially more luminous, but that 
were subsequently decimated by disk (or core) shocks and subsequent tidal detachment of stars. This hypothesis is 
supported by the observation that the faint cluster E3 (My = —2.8, R = 7.6 kpc) is unusually rich in binary stars 
(van den Bergh 1980). 
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3. CONCLUSIONS 

The following conclusions may be drawn from the presently available data on Galactic globular 
clusters: 

• Collapsed cores are rare among the most luminous Galactic globulars. This might suggest 
that the the collapse time scale for the most massive globular clusters is longer than the age 
of the Galaxy. 

• Most of the globular clusters that do not have collapsed cores show a strong correlation 
between luminosity and central concentration. This correlation is in the sense that the most 
luminous clusters have the strongest central concentration of light. 

• Metal-poor clusters with [Fe/H] < —1.0 occur at all Galactocentric distances, whereas all 
but three of the metal-rich clusters are situated at R < 10 kpc. 

• These three metal-rich ([Fe/H] > —1.0) clusters at R > 10 kpc may originally have formed 
in, or in association with, dwarf spheroidal companions to the Galaxy. 

• For metal-rich clusters with [Fe/H] > — 1.0 there is no evidence for the existence of a radial 
metallicity gradient. The data hint at, but are not numerous enough to establish, the reality 
of a radial metallicity gradient among metal-poor clusters with [Fe/H] < —1.0. 

• Both the outer halo globular clusters with R > 80 kpc, and those associated with the Sagit- 
tarius dwarf spheroidal, appear to have bi-modal luminosity functions. It is difficult to see 
how such bi-modal luminosity distributions could have arisen from clusters that initially had 
a single power law mass spectrum. 

The luminosity My, the metallicity [Fe/H], the Galactocentric distance R, and the concentra- 
tion parameter c, are presently known for the great majority of Galactic globular clusters. Never- 
theless, the relatively small total number of Galactic globular clusters places significant statistical 
restrictions on the degree of certainty with which correlations between many of these parameters 
can be established. In particular the Galactic globular cluster population is too small to establish if 
there exists a radial metallicity gradient among the metal-poor ([Fe/H] < —1.0) halo cluster pop- 
ulation. By the same token it is not possible to say if the metal-rich bulge/disk globular clusters 
with [Fe/H] > — 1.0 have a radial composition gradient. 

It is a pleasure to thank Russell Redman for his kind help with Figure 1. I also thank Scott 
Tremaine and the referee for helpful comments. 
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Table 1. Galactic globular cluster parameters 
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Table 1 — Continued 
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Table 1 — Continued 
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N6362 


-0.95 


1.10 


-7.06 


0.72 


Lil 1 


+0.22 


2.30 


-7.63 


0.41 


N6380 


-0.50 


1.55 


-7.46 


0.51 


Ter 1 


-0.35 


2.50 


-3.30 


0.26 


Ton 2 


-0.50 


1.30 


-6.14 


0.15 


N6388 


-0.60 


1.70 


-9.82 


0.64 


N6402 


-1.39 


1.60 


-9.02 


0.59 


N6401 


-1.12 


1.69 


-7.62 


-0.10 


N6397 


-1.95 


2.50 


-6.63 


0.78 


Pal 6 


-0.10 


1.10 


-7.37 


-0.10 


Dj 1 




1.50 


-6.40 


0.00 


Ter 5 


-0.28 


1.74 


-7.91 


-0.15 


N6440 


-0.34 


1.70 


-8.75 


0.11 


N6441 


-0.53 


1.85 


-9.47 


0.54 


Ter 6 


-0.50 


2.50 


-7.87 


0.20 


N6453 


-1.53 


2.50 


-7.05 


0.52 


UKS 1 


-0.50 


2.10 


-6.88 


-0.10 


N6496 


-0.64 


0.70 


-7.23 


0.63 


Ter 9 


-1.00 


2.50 


-3.93 


-0.22 


Dj 2 


-0.50 


1.50 


-6.98 


0.15 


N6517 


-1.37 


1.82 


-8.28 


0.63 


Ter 10 


-0.70 




-6.31 


0.38 
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Table 1 — Continued 



Cluster 


\Fe/H] 

L / J 


c 


My 


log-R 


N6522 


-1.44 


2.50 


-7.67 


-0.22 


N6535 


-1.80 


1.30 


-4.73 


0.59 


N6528 


-0.17 


2.29 


-6.93 


0.11 


N6539 


-0.66 


1.60 


-8.30 


0.49 


N6540 


-1.20 


2.50 


-5.38 


0.64 


N6544 


-1.56 


1.63 


-6.56 


0.73 


N6541 


-1.83 


2.00 


-8.37 


0.34 


N6553 


-0.34 


1.17 


-7.99 


0.40 


N6558 


-1.44 


2.50 


-6.46 


0.00 


11276 


-0.73 


1.29 


-6.67 


0.57 


Ter 12 


-0.50 


0.57 


-4.14 


0.53 


N6569 


-0.86 


1.27 


-7.88 


0.08 


N6584 


-1.49 


1.20 


-7.68 


0.85 


N6624 


-0.42 


2.50 


-7.50 


0.79 


N6626 


-1.45 


1.67 


-8.33 


0.41 


N6638 


-0.99 


1.40 


-6.83 


0.20 


N6637 


-0.71 


1.39 


-7.52 


0.20 


N6642 


-1.35 


1.99 


-6.57 


0.20 


N6652 


-0.96 


1.80 


-6.57 


0.38 


N6656 


-1.64 


1.31 


-8.50 


0.69 


Pal 8 


-0.48 


1.53 


-5.52 


0.75 


N6681 


-1.51 


2.50 


-7.11 


0.32 


N6712 


-1.01 


0.90 


-7.50 


0.54 


N6715 


-1.59 


1.84 


-10.01 


1.29 


N6717 


-1.29 


2.07 


-5.67 


0.36 


N6723 


-1.12 


1.05 


-7.86 


0.41 


N6749 


-1.60 


0.83 


-6.70 


0.70 


N6752 


-1.56 


2.50 


-7.73 


0.72 


N6760 


-0.52 


1.59 


-7.86 


0.68 


N6779 


-1.94 


1.37 


-7.38 


0.99 


Ter 7 


-0.58 


1.08 


-5.05 


1.20 


Pal 10 


-0.10 


0.58 


-5.79 


0.81 


Arp 2 


-1.76 


0.90 


-5.29 


1.33 
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Table 1 — Continued 



Cluster 


[Fe/H] 


c 


My 




N6809 


-1.81 


0.76 


-7.55 


0.58 


Ter 8 


-2.00 


0.60 


-5.05 


1.28 


Pal 11 


-0.39 


0.69 


-6.81 


0.89 


N6838 


-0.73 


1.15 


-5.56 


0.83 


N6864 


-1.32 


1.88 


-8.35 


1.11 


N6934 


-1.54 


1.53 


-7.65 


1.16 


N6981 


-1.40 


1.23 


-7.04 


1.11 


N7006 


-1.63 


1.42 


-7.68 


1.59 


N7078 


-2.25 


2.50 


-9.17 


1.02 


N7089 


-1.62 


1.80 


-9.02 


1.02 


N7099 


-2.12 


2.50 


-7.43 


0.85 


Pal 12 


-0.94 


1.94 


-4.48 


1.20 


Pal 13 


-1.65 


0.66 


-3.51 


1.44 


N7492 


-1.51 


1.00 


-5.77 


1.40 
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Table 2. Radial dependence of metallicity 





n 




n 




logi? 


-2.40 < [Fe/H] 


< -1.60 


-1.60 < [Fe/H] 


< -1.00 


2.00 to 2.19 


1 


1 


1 


1 


1.80 to 1.99 


2 


2 


2 


2 


1.60 to 1.79 








1 


1 


1.40 to 1.59 


6 


6 


3 


3 


1.20 to 1.39 


11 


6 


5 


3 


1.00 to 1.19 


5 


5 


6 


6 


0.80 to 0.99 


8 


8 


6 


6 


0.60 to 0.79 


3 


3 


9 


9 


0.40 to 0.59 


5 


5 


7 


7 


0.20 to 0.39 


4 


4 


6 


6 


0.00 to 0.19 


2 


2 


2 


2 


-0.20 to -0.01 








2 


2 


-0.40 to -0.21 








2 


2 



1 Numbers in fourth and sixth columns of the table exclude clusters as- 
sociated with the Sagittarius system. 
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Table 3. Magnitude dependence of cluster central concentration 



My 




1.50 < c < 2.00 


c > 2.40 


-10.00 to ■ 


-10.99 


1 





-9.00 


-9.99 


7 


1 


-8.00 


-8.99 


16 


1 


-7.00 


-7.99 


15 


12 


-6.00 


-6.99 


5 


5 


-5.00 


-5.99 


1 


2 


-4.00 


-4.99 


1 





-3.00 


-3.99 





2 


-2.00 


-2.99 


1 






